Weak nonmild solutions to some SPDEs* 



Daniel Conus 
University of Utah 



Davar Khoshnevisan 
University of Utah 



April 15, 2010 



Abstract 



We study the nonlinear stochastic heat equation driven by space- 
time white noise in the case that the initial datum uq is a (possibly 
signed) measure. In this case, one cannot obtain a mild random-field 
solution in the usual sense. We prove instead that it is possible to 
establish the existence and uniqueness of a weak solution with values 
in a suitable function space. Our approach is based on a construction 
of a generalized definition of a stochastic convolution via Young-type 
inequalities. 

Keywords: Stochastic PDEs, weak solutions, measure-valued initial 
conditions. 
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1 Introduction 

Let us consider the nonlinear stochastic heat equation 



— u t (x) = (Cu t )(x) + a(u t (x))W(t ,x) (t > 0, x G R), (1.1) 



where: (i) C is the generator of a real- valued Levy process {X t }t>o with 
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Levy exponent \E r , normalized so that Ee 4< ^* 



e for every £ £ R and 



t > 0; (ii) a : R — > R is Lipschitz continuous with Lipschitz constant Lip CT ; 
(hi) W is space-time white noise; and (iv) The initial datum uq is a signed 
Borel measure on R. 

Equation (1.1) arises in many different contexts; three notable examples 
are Bertini and Cancrini [1], Gyongy and Nualart [15], and Carmona and 
Molchanov [6]. 

In the case that uq : R — > R+ is a bounded measurable function, the 
theory of Dalang [10] shows that there exists a unique random-field mild 
solution {ut(x)}t>o,xeR provided that 

1 f°° df 

W): =^L p + 2Re m <0 ° ^ some, hence all, /3 > 0. (1.2) 

In general, Dalang's Condition (1.2) cannot be improved upon [10, 19]. 

Dalang's condition (1.2) implies also that the Levy process X has tran- 
sition functions pt(x) [14, Lemma 8.1]; i.e., for all measurable / : R — > R + , 

/oo 
p t (z)f(z)dz forallt>0. (1.3) 
-oo 

A mild solution in this setting is a random field {ut(x)}t>o,x<=R that satisfies 

u t {x) = (P t u )(x) + / Pt- S {y - x)a(u s {y))W(dsdy) a.s., (1.4) 
J[0,i]xR 

for all t > and x G R, where {Pt}t>o denotes the semigroup associated 
to the process X, and the stochastic integral is understood as a Walsh 
martingale-measure stochastic integral [20]. Notice that (1.4) can be rewrit- 
ten in the following form: For all (t ,x) G R + x R, 

ut(x) = (P t u )(x) + (p* (a o u)W^j (x) a.s., (1.5) 

where "*" denotes "stochastic convolution"; see (2.1) below and pt(x) = 
Pt(—x) for all x G R. 

In the case that uq is not a bounded and measurable function, but instead 
a (possibly signed) Borel measure on R, the solution u cannot be defined as 
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a random field, but has to be considered as a process of function-space type. 
As a consequence, the stochastic convolution in (1.4) is not well defined in 
the sense of Walsh. Section 2 below is devoted to extending the definition 
of the stochastic convolution of a process T with respect to ZW in the case 
that Z is in a suitable Banach space Bp of random processes. The key step 
of this extension involves developing a kind of "stochastic Young inequality" 
(Proposition 2.1). Such an inequality appeared earlier in [9], in a different 
context, in order to obtain intermittency properties for equation (1.1) in the 
case that uq is a lower semicontinuous bounded function of compact support. 

In Section 3 we establish the existence and uniqueness of a weak solution 
to (1.1). Namely, we prove that Dalang's condition (1.2) implies that if 
uq = ii is a (possibly signed) Borel measure on R that satisfies a suitable 
integrability condition (3.5), then there exists a unique u G Bp such that u 
almost surely satisfies (1.5) for almost every t > and x € R. This solution 
is not a random field; but it is a function-space-valued solution. 

In Section 4 we prove that our condition for existence and uniqueness is 
unimprovable. And in Section 5 we mention briefly examples of initial data 
uq that lead to the existence and uniqueness of a weak solution to (1.1), 
together with further remarks that explain what happens if we study the 
1-D stochastic wave equation in place of (1.1). 

2 Generalized stochastic convolutions 

Let T : (0,oo) x R — > R be measurable, and Z := {Z t (x)} t >o,x£Ti be a 
predictable random field in the sense of Walsh [20, p. 292]. Let us define the 
stochastic convolution V * ZW of the process T with the noise ZW as the 
predictable random field 

(T*ZW) t (x):= f T t _ s (x-y)Z s (y)W(dsdy). (2.1) 
J[0,i]xR 

The preceding is defined as a stochastic integral with respect to the mar- 
tingale measure ZW in the sense of Walsh [20, Theorem 2.5], and is well 
defined in the sense of Walsh [20, Chapter 2] provided that the following 



3 



condition holds for all t > and x G R: 

2 rt 



2 ./o 



(T*m) ( (i) = / ds I dy [T t _ s {x-y)} 2 \\Z s {y)\\l < oo. (2.2) 



Let W 2 denote the collection of all predictable random fields Z that 
satisfy the following: Z t (x) G L 2 (P) and 

ft roo 

ds dy [F T _ s (x - y)} 2 \\Z s (y)g < oo, (2.3) 

JO J-oo 

for all < t < t and x G R. 

We may think of the elements of W 2 as Walsh-integrable random fields. 
And because (2.3) implies (2.2), the preceding discussion tells us that the 
stochastic convolution V * ZW is a well-defined predictable random field for 
every Z G W 2 . 

Our present goal is to extend the definition of the stochastic convolution 
of Z so that it includes more general random processes Z. Other extensions 
of this stochastic convolutions have been developed for other purposes as 
well [8,10,11,17]. 

Let us choose and fix a real number k G [2 , oo), and define L k to be the 
collection of all predictable random fields {Zt(x)}t>o,xe~R, such that Zt(x) G 
L fe (P) for all t > and x G R. Let M(R) be the space of (T-finite Borel 
measures on R. For every j3 > and 77 G M(R) we can define a norm on 
L k as follows: For every Z G L k , 

/ f-oo roo \ 1/2 

A/^Z) := / e~*dt sup / V (dx) \\Z t (x - z)\\ 2 k ) . (2.4) 

Here and throughout, we use implicitly the following observation: If 
Z, Z' G L k satisfy M k v (Z - Z') = 0, then Z and Z' are modifications of one 
another. There is an obvious converse as well: If Z and Z' are modifications 
of one another, then Np v (Z — Z') = 0. We omit the elementary proof. 

Our next proposition is a "stochastic Young's inequality," and plays a 
key role in our extension of Walsh-type stochastic convolutions. But first 
let us introduce some notation that will be used here and throughout. 
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Throughout this paper, Zk denotes the optimal constant in Burkholder's 
L fc (P)-inequality for continuous square-integrable martingales; its precise 
value involves zeroes of Hermite polynomials, and has been computed by 
Davis [13]. By the ltd isometry, Z2 = 1. Carlen and Kree [5, Appendix] 
have shown that Zk < 2a/& for all k > 2, and moreover Zk = {2 + o(l))Vk 
as k — > oo. 

Proposition 2.1 (A stochastic Young's inequality). For every k € [2,oo), 
Z eW 2 nL k ,r]e Af (R), and /3 > 0, 

/ /-oo \ 1/2 

AA^(r*^)<z fe ^ e-^||r t ||| 2(R) dtj -AA^(Z), (2.5) 

where is the optimal constant in Burkholder's inequality. 

Remark 2.2. We emphasize that W 2 n L 2 = W 2 . □ 

Before we prove Proposition 2.1 let us first describe how it can be used 
to extend stochastic convolutions. Proposition 2.1 will be proved after that 
extension is described. 

Let Bp denote the completion of W 2 n L k under the norm A/"^. 1 It 
follows then that is a Banach space of predictable processes [identified 
up to evanescence]. 

Proposition 2.1 immediately implies that if 

POO 

T(/3):= ^ e-^||r t ||| 2(R) dt<oo, (2.6) 

then Z i— > T * ZW has a unique extension to all Z G B k , and the resulting 
extension — written still as Z \— > T * ZW — defines a bounded linear operator 
from B k into itself. And the operator norm is at most the square root of the 
Dalang integral T(/3). [In the case that Tj(x) denotes the transition density 
of a Levy process with Levy exponent ^, Plancherel's theorem implies that 
T(/3) is the same Dalang integral as in (1.2); see (3.1) below as well.] 

lr rhe latter is of course a norm on equivalence classes of modifications of random fields 
and not on random fields themselves. But we abuse notation as it is standard. 
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From now on, we deal solely with this extension of the stochastic convo- 
lution. However, we point out the there is a great deal of variability in this 
extension, as the parameters /3 > 0, k G [2 , oo), and r\ G M(R) can take on 
many different values. 

Let us conclude this section by establishing our stochastic Young's in- 
equality. 

The proof of Proposition 2.1 relies on an elementary estimate for Walsh- 
type stochastic integrals. 

Lemma 2.3. For all real numbers t > 0, x G R, and k G [2,oo) ; and for 
every Z G W 2 n L k , 

\(r*ZW) t (x)\\ k <4[j Q ds J Jy [r t . s (x-y)] 2 \\Z s (y)\\l) • (2-7) 
Proof. Condition (2.3) implies that if < t < r, then 



(r * ZW\ T {x) := / IV s (x - y)Z s (y) W(ds dy) (2. 
'[0,t]xR 



is well defined and in L 2 (P). Moreover, 

(r*w)i, T (i) 



t roo \ I/ 2 

ds / dy [T T ^{x-y)} 2 \\Z s {y)\\ 2 ) . (2.9) 

J-oo J 



Walsh's theory of martingale measures [20, Theorem 2.5] tells us that the 
stochastic process (0,r] 3 t >->■ (r * ZW)t, T {x) is a continuous L 2 (P)- 
martingale. Therefore, Davis's refinement [13] of Burkholder's inequality 
[2, 3, 16] implies that 



(T*ZW\ T {x) 



<4 



/oo 
dy [T T _ s (x - y)} 2 [Z s (y)Y 
-oo 



k/2 



k/2 



(2.10) 



And it follows from Minkowski's inequality that 



(T*ZW) t , T (x) k k <z k k j\s I dy[T T ^(x-y)Y\\Z s {y)f k . (2.11) 
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The lemma follows from this upon setting r := t. 



□ 



Proof of Proposition 2.1. The original construction of Walsh implies that 
|| (r * -ZW)t(x)||fc defines a Borel- measurable function of (t , x) £ (0 , oo) x R. 
Indeed, it suffices to verify this measurability assertion in the case that Z is 
a simple function in the sense of Walsh [20, p. 292], in which case the said 
measurability follows from a direct computation. 

We may apply Lemma 2.3 with x — z in place of the variable x, and then 
integrate [drj] to obtain 



/oo 
r)(dx) (T*ZW) t (x-z) 
-CO 

/OO f't f'OO 

rj(dx) ds dy [r t _ s 0r - z - y)f\\Z s {y)\\l 
-oo JO J —oo 

/OO f't f'OO 

n(dx) ds dy [T t . s {y)?\\Zs(x - z - y)f k 
-oo JO J —oo 

f't f'OO 

< z\ ds ||rV s ||2 2(R)S up / r](dx) \\Z s (x - v)\\ 2 k . 



(2.12) 



veR J — oo 



Or equivalently, 

/oo 
rj(dx) (T*ZW) t (a 
-oo 

ft f'OO 

<z\\ ds ||r t _ s ||| 2(R) sup / n(dx) \\Z s (x - z)\\ 2 k . 

JO zeRJ-oo 



(x - z, 



(2.13) 



Multiply both sides by exp(— /?£), integrate [dt] and use Laplace transforms 
properties for convolutions to obtain the result. □ 

Proposition 2.4. Suppose a : R — > R is Lipschitz continuous and Z, Z* € 
.Bg for some k £ [2 , oo), f3 > 0, and r] G M(R). T/ien, 

o Z - a o Z*) < Li P(T • AT| ;7? (Z - Z*). (2.14) 

Proof. If Z, Z* G n -L fc , then this is immediate. In the general case 
we proceed by approximation: Let Z 1 , Z 2 , . . . , Z 1 '*, Z 2 '* , ... be in W 2 n L fc 
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such that Z n -»• Z and Z n '* -»• Z* in Bj >tJ , as n -> oo. By going to a 
subsequence, if necessary, we can [and will!] also assume that 

Af^ v {Z n - Z n+1 ) + Af^ v {Z n '* - Z n+1 '*) < 2- n for all n > 1. (2.15) 

It follows also that for all n > 1, 

Ml^a oZ n -ao Z n+l ) +Nl v {a o Z n '* - a o Z n+1 >*) < Li P(7 • 2-". (2.16) 

Of course, this implies immediately that a o Z n and a o Z n> * converge in 
. It suffices to prove that the mentioned limits are respectively a o Z 
and a o Z* . But (2.16) implies that 



. a{Z?(x))-a{Zr\x))-or 
. a{Zr{x))-a{Zr^{x)). 
Because Yl n =i su Pzgr( " ' ) — su PzeR Yl^Lii ' ' ' ) i n (2-17), it follows readily 



that for almost every pair (t ,x) G R + x R: 

• lim n _ i , 00 Z™(x) = Z t {x) almost surely; 

• limn-s.oo Z™'*(x) = Zf(x) almost surely; 

• lim^oo <r(Zt(x)) = a(Z t (x)) almost surely; and 

• lim„^ 00 a(Zf'*(x)) = a{Z^{x)) almost surely. 

[Note the order of the quantifiers!] We showed earlier that lim ra _ >00 a o Z n 
and limn^oo a o Z n '* exist in . The preceding shows that those limits 
are respectively a o Z and a o Z* . This completes the proof. □ 




(2.17) 



where A™(x) stands for any one of the following four quantities: 
. Z?{x)-Z? +l {x)- 

. zr(x)-Z? +1 '*(x); 
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3 Existence and Uniqueness 



This section is devoted to the statement and proof of the existence and 
uniqueness of a weak solution to (1.1). We will make use of the generalized 
stochastic convolution developed in Section 2. 

Before we proceed further, let us observe that from now on T t (x) of the 
previous section is chosen to be equal to the modified transition functions 
Pt{x), in which case Dalang's integral can be computed from PlancherePs 
formula as follows: 

poo 1 poo J( 

In particular, the T of (2.6) and that of (1.2) are equal in the present setting. 

Next we identify our notion of "solution" to (1.1) in the case that uq = fx 
is a measure. 

Suppose first that uq is a nice function and (1.1) has a mild solution u 
with initial datum no- Then for alH > and x 6 R, 

P [u t (x) = (P t u )(x) +(p*(ao u)w) (*)} = 1. (3.2) 

Consequently, Fubini's theorem tells us that every mild solution u to (1.1) 
with initial function fj, := no is a weak solution in the sense that the following 
holds with probability one [note the order of the quantifiers!]: 

ut(x) = (Ptv)(x) + (p* (job)?) (i) for a.e. (t , x) G R+ x R. (3.3) 

It is easy to see that the preceding agrees with Walsh's definition of a weak 
solution [20, p. 309]. 

Now we consider (1.1) in the case that uq = u is a possibly-signed Borel 
measure on R. 

Let us suppose that Dalang's condition (2.6) holds, and consider an 
arbitrary u G . Since a is Lipschitz continuous, it follows that a o u G 
Bp . Therefore, we can conclude that the stochastic convolution p*(aou)W 
is a well-defined mathematical object, as was shown in the previous section. 
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Consequently, we can try to find a solution u to (1.1) with uq = n by seeking 



to find u G Bp such that 

n = P,fi + p* (aou)W, (3.4) 

where the equality is understood as equality of elements of Bp . Of course, 
we implicitly are assuming that P #/ u G #g as well. That condition is clearly 
satisfied if 



/"OO /"OO 

/ e-^dssup/ 7?(dx)|(P s/ u)(x-z)| 2 < oo. (3.5) 

JO zeRJ-oo 

Then, u is a solution of function-space type to (1.1) with no = «. But it 
has more structure than that. Indeed, suppose that: (i) (3.5) holds; and (ii) 
There exists u G Bp that satisfies (3.4). Then the preceding discussion 
shows also that u is a weak solution to (1.1) in the sense of Walsh [20, p. 
309]. And it would be hopeless to try to prove that such a u is a mild 
solution, as there is no natural way to define ut(x) for all t > and x G R. 

Throughout the remainder of this section we choose n G M(R). In the 
case that <r(0) / 0, then we assume additionally that r] is a finite measure. 

Theorem 3.1. Consider (1.1) subject to uq = fi, where u is a signed mea- 
sure that satisfies (3.5). // (1.2) holds and 

then there exists a solution u G Bp that satisfies (3.4) . Moreover, u is 
unique in Bp^; i.e., if there exists another weak solution v that is in Bp^ 
for some k > 2, then v is a modification of u. 

Proof. First, we argue that we can always choose j3 such that (3.6) holds. 

Indeed, Condition (3.5) implies that Mp v (P 9 uq) < oo for all (3 > and 
k G [2 , oo). Also, because of Dalang's condition (1.2), and by the monotone 
convergence theorem, lim Q _ >00 T(a) = 0. Therefore, we can combine these 
two observations to deduce that (3.6) holds for all (3 large, where 1/0 := oo. 
Throughout the remainder of the proof, we hold fixed a /3 that satisfies (3.6). 
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Set := 0, and iteratively define 

> +1 ) := P.^ + p* ([a o U W]W) . (3.7) 



These u^ n+1 ^s are all well denned elements of Bp . In fact, it follows from 
Proposition 2.1 that for all n > 0, 

Ml v (> +1 )) < Ml {P.n) + z k ^W)^l (a o „W) . (3.8) 

And because |<r(z)| < |<r(0)| + Lip CT |z| for all zeR, 



< Ml (P.m) + **>/TG9) [k(0)| • A/^l) + Lip CT • Ml («<»>) 



(3.9) 



where := 1 for all £ > and x G R. In particular, G 5^ for all 

^ > 0, by induction. This is clear if <r(0) = 0; and if <r(0) ^ 0, then it is also 
true because Ml(l) = y^R)//? < oo, thanks to the finiteness assumption 
on rj [for the case <r(0) / 0]. Moreover, (3.6) and induction together show 
more; namely, that sup„> Mp „(u^ ) < oo. 

A similar computation, this time using also Proposition 2.4, shows that 
for all n > 1, 

Ml (> +1 ) - «(»>) < z fc Li P(T yrM • Ml («W - n^- 1 )) . (3.10) 

And (3.6) implies that £^JV$ >TJ (u (n+1) -^ (n) ) < °o, therefore, {u^}%> =0 
is a Cauchy sequence in Bp v - Let u := hin^^oo , where the limit takes 
place in . According to Proposition 2.1, 



Ml (p * uWw - p * u#) < z k y/T(P) ■ Ml (k W - « 



(3.11) 



as n — > oo. 



It follows readily from these remarks that Mp(u — P,/j, + p*(ao u)W) = 0. 
That is, u satisfies (3.3) for almost all (t,x) G R + x R; see also (3.4). This 



11 



proves part (1) of the theorem. 

In order to prove the second part, let us suppose that there exists another 
"weak solution" v G B% . Then, 5 := u - v G B% and 

$ = P * (\? ° u]W^ — p* ([a o = p * ^[cr oii-(To v]W^j . (3.12) 

[The second identity follows from the very construction of our stochastic 
convolution, using the fact that Z p* ZW is a bounded linear map from 
Bp v to itself.] Propositions 2.1 and 2.4 together imply the following: 

Nl v {S) < Z k ^W)- Kr, (<T°U-aav) 

Thanks to (3.6), .A/J (u — u) = A/"^ (<5) = 0. This readily implies that u and 
v are modifications of one another, as well. □ 



4 On Condition (3.6) 

Let us consider the measure ry m G M(R) defined by 

Vm (dx) = e~^ m dx, (4.1) 

where m > is fixed. If <r(0) = 0, then we may take m := oo, whence 
rj(dx) = dx. 

Theorem 4.1. Suppose (1.1) /ias a solution u G n m >o-Bjg wrai/j uq = fx 
for a nonvoid signed Borel measure fi on R iwi/t |/u|(R) < oo. Suppose also 
that L CT := inf^ € R |cr(£)/z| > 0. Then, /3 satisfies T(/3) < L~ 2 . 

Proof. Let A^^ be the norm defined by 

(/•oo /<oo \ 1/2 

e-^dt y e-W/ m dx \\Z t (x)\\lj . (4.2) 

Notice that Mp is similar to Np^m > but is missing a supremum on Z in 
the space variable; cf. (2.4). Moreover, Mp(u) < ftfp Vm (u) < oo. Note 
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that if H, Z e B 2 R „ with one of them — say H — random and the other 
one deterministic, then we have [Mp(H + G)} 2 = [Mp(H)] 2 + [Mp{G)\ 2 . 
This is a direct computation if H, G € W 2 ; the general case follows from 
approximation (we omit the details because the method appears already 
during the course of the proof of Proposition 2.4). It follows that 

[Mp(u)] 2 = [Mp(P.fj)] 2 + [m p (p*([a o u}W))] 2 . (4.3) 
The method of proof of Proposition 2.4, together with the simple bound, 

e -\x\/m > e -\x-y\/m . e -|tf|/m > s h OW S also that 

Mp(p*([tTou]W)^>L a Mp(p*uWy (4.4) 

But 

/ roo \ 1/2 

Mp(p*ZW)=U e-^\\p t \\ 2 LHn) dt) -Mp(Z). (4.5) 
[Again one proves this first for nice Z's and then take limits.] Therefore, 

Mp (p * {[a o u ]wj) > UVW)' M p {u). (4.6) 

Combine this with (4.3) to find that 

[Mfs(u)} 2 > [Mfs(P.fi)} 2 + L 2 a T(P) [Mp{u)} 2 . (4.7) 

Now suppose, to the contrary, that T(/3) > L^ 2 . Then, it follows that 
Mp(Pmfi) = regardless of the value of m; i.e., for all m > 0, 

/*oo />oo 

/ e-^dt e-W m dx \{Pt^){x)\ 2 = 0. (4.8) 

JO J-oo 

Let m | oo and apply the monotone convergence theorem, and then the 
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Plancherel theorem, in order to deduce that 

poo 

= J e-^||P t H|| 2(R) d* 



f'OO f'OO 

/ e-^dt dee^ Rc *«)|/i(0| 2 (4.9) 

JO J-oo 



2vr . , 

i r iacoi 2 , de 



Since ^ is never infinite, the preceding implies that fx = 0, which is a 
contradiction. It follows that T(/5) < L~ 2 . □ 

Theorem 4.1 implies also that Condition 3.6 is sharp: Consider the case 
that Lip CT = L CT . [This is the case, for instance, for the parabolic Anderson 
problem where a(x) oc x, or when j3 has sharp linear growth.] Then in this 
case Theorem 3.1 and Theorem 4.1 together imply that (3.6) is a necessary 
and sufficient condition for the existence of a weak solution to (1.1) that has 
values in C\ m> iB k B . 

5 Examples and Remarks. 

Example 5.1 (A parabolic Anderson model). Let cr(x) = Xx. In that case, 
the solution u corresponds to the conditional expected density at time t > 
of a branching Levy process starting with distribution uq , given white- noise 
random branching. The case that a(0) =0 and uq is a function with compact 
support is studied in [9], in which intermittency properties are derived. Here, 
uq can be a compactly supported measure (not necessarily a function). If 
we let i] denote the one-dimensional Lebesgue measure, then (3.5) becomes 



f 

j — < 



l''«>' 2 d f <oo. (5.1) 



/3 + 2Re*(£) 



For instance, if [i = 8q, then condition (5.1) is precisely Dalang's condition 
(1.2), and (1.1) admits a weak solution. In this way we can now define 
properly the solution of the parabolic Anderson model with u$ = 5q, which 
was studied in Bertini and Cancrini [1] . □ 
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Remark 5.2 (A nonlinear stochastic wave equation). It is possible to apply 
similar techniques to the study of the following nonlinear stochastic wave 
equation driven by the Laplacian: 



j^ut(x) = K 2 (Au t )(x) + a(u t (x))W(t,x) (t>0,z€R). (5.2) 



If the initial conditions uq and vq are nice functions, then the solution to 
(5.2) can be written as 



i[0,t]xR 

where V is the fundamental solution of the 1-dimensional wave equation, 
namely 



and T' t denotes the weak spatial derivative of iy Then, the existence and 
uniqueness of a weak solution to (5.2) in the case that uq and vq are (possibly 
signed) Borel measures on R can be established using the techniques of this 
paper, since the definition of the generalized stochastic convolution applies 
as well with the 1-D wave propagator T above. The conditions on the initial 
conditions have to be adapted to insure that T' t * uq and T t * vq both are 
in Bp , but are similar to (3.5). The details on this are left to the reader, 
as the stochastic wave equation in dimension one has been widely studied 
already [4,7, 12,18-21]. 
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